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We consider the Color Glass Condensate (CGC) at finite temperature, which would be relevant 
to the initial condition of relativistic heavy ion collisions and the energy loss of energetic partons in 
the quark-gluon plasma. In the weak source approximation, we derive the thermal BFKL equation. 
We find that the thermal effect shows up as a Bose enhancement of the soft gluon emission if the 
temperature is high enough to match the light-cone energy of the soft gluons. This suggests that 
the saturation regime could be reached sooner than in the vacuum. 
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I. INTRODUCTION 

The Relativistic Heavy Ion Collider (RHIC) at BNL 
has revealed various interesting natures of high energy 
QCD matter which would be the quark-gluon plasma 
(QGP). One of the findings is that the expanding mat- 
ter can be described as a perfect fluid without dissipa- 
tions P, 0, 0, 0|- Even surprising would be the early 
thcrmalization. Hydrodynamic analyses suggest that 
the thermalization takes place in less than 1 fm/c af- 
ter the collision 2]. In the future Large Hadron Col- 
lider (LHC) experiments at CERN, where higher en- 
ergy collisions are possible, we expect an even earlier 
thcrmalization, that could be almost instantaneous af- 
ter the collision moment @. It is important issues 
to elucidate the mechanism of the early thermalization 
[|, 0, H i, US El El [3 , as well as to determine the ini- 
tial condition to be used as an input of hydrodynamical 
simulations @, EI EE E3, EE El ■ Another discovery 
to be mentioned would be the jet quenching. QGP has an 
enormous stopping power for energetic partons, for which 
full theoretical understanding has not been achieved. 

On the other hand, small x physics in high energy 
hadron scatterings is one of the fascinating areas in 
QCD physics. The gluon distribution in the small x re- 
gion, which is responsible for the high energy hadronic 
scattering amplitude, is determined by the BFKL equa- 
tion [2(| [U [22j]. It is known, however, that the 
BFKL equation predicts too many gluons as we go to 
smaller x, resulting in unitarity violation of the scatter- 
ing cross-section. We need some saturation mechanism 
to attenuate the growth of the small x gluon number 
[EH Hllllil- The Color Glass Condensate (CGC) is 
one of the formalisms which successfully implements the 
saturation [27|, S [H [H M, M, H, 0, H M 13 • The 
idea of CGC is simple and clear. We assume the pres- 
ence of recoilless large x partons which travel at the speed 
of light and stay on the light-cone. These large x par- 
tons act as static color sources, which radiate smaller x 
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partons. And the emitted gluons radiate further smaller 
x partons. At some point, however, the population be- 
comes so large that the recombination between gluons 
begins to be effective, resulting in the saturation in the 
small x region. The state in the saturated regime, where 
the gluon field can be treated as classical, is called the 
Color Glass Condensate. Mathematically, the successive 
emission of smaller and smaller x gluons is formulated 
as a Renormalization Group (RG) equation of the effec- 
tive action. Within the formalism, the BFKL equation 
can be reproduced in the weak field regime. In the satu- 
ration regime, the evolution equation is replaced by the 
JIMWLK equation, which is a non- linear extension of the 
BFKL equation. 

What we intend to do in this paper is to extend the 
CGC formalism to finite temperature and to derive the 
thermal BFKL equation. There are two reasons for do- 
ing that. (1) As we have mentioned, the matter reaches 
thermalization just after the collision, through some very 
complicated processes. Actually, the mechanism of the 
early thcrmalization is not fully understood so far. It 
may be that the mechanism should be verified within 
the CGC formalism itself [H, El- For our purpose, 
let us assume the early thermalization as an empirical 
fact. On the other hand, in the Bjorken picture of the 
ultra-relativistic nuclear collision [3a ]. the large x par- 
tons including valence partons can go through the other 
nucleus without damage, and remain almost on the light- 
cone. Then, just after the collision, the large x partons 
would radiate smaller x partons in the medium at finite 
temperature. This gluon distribution, as is determined 
by the thermal BFKL or JIMWLK equation, would give 
us the initial condition for the hydrodynamical expan- 
sion that follows. It may be that CGC is reconstructed 
at that moment, which is not necessarily identical with 
that being present before the collision [39j. In fact, in 
Refs. [10, EH, EH , it was shown that CGC can become an 
appropriate initial condition. 

(2) The other motivation is related to the jet quenching 

[43l [44j . As is discussed in Ref. (44| , the energetic par- 
ton going through QGP can be regarded as a collision of 
the parton and the QGP with the large center-of-mass 
energy. Just as in the high energy hadronic scattering 
where the hadrons exhibit the saturation, the QGP would 
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have saturated gluons as a consequence of radiation by 
the thermal quarks and gluons. Those extra emitted glu- 
ons could explain the observed large stopping ability of 
QGP. We notice that this argument is in a frame where 
the thermal bath of QGP is moving at high speed. In the 
rest frame of QGP, we have another picture, which should 
be equivalent to the physics in the QGP moving frame. 
Here the energetic parton radiates soft gluons inside the 
medium at finite temperature that is at rest. Those ra- 
diated soft gluons, the distribution of which could be de- 
termined by the thermal BFKL equation derived in the 
CGC formalism, would be responsible for the parton's 
energy loss in QGP. 

Motivated by these reasons, we will extend the CGC 
formalism to finite temperature to derive the thermal 
BFKL equation. 

The organization of this paper is as follows. In Sec. Hi] 
we present the gluon propagator in the light-front field 
theory at finite temperature. In Sec. IIII1 we discuss the 
CGC formalism at finite temperature. In Sec. IIV1 we 
derive the thermal BFKL equation. Summary is given in 
Sec.|V] 

Throughout this paper, we use the same convention of 
the light-cone variable and the metric and so on as in 
Ref. [3|. 



II. LIGHT-FRONT FIELD THEORY AT FINITE 
TEMPERATURE 

Let us begin with basic concepts of (conventional) 
thermal field theory [45], [46|. All the thermodynamic 
quantities can be derived from the partition function 
Z = Trexp(— (3H) with the Hamiltonian H and the 
inverse temperature (5, which can be rewritten in the 
path integral form. There are two formalisms avail- 
able: The imaginary-time formalism and real-time for- 
malism. In the imaginary-time formalism, the (anti- 
periodic boundary condition with respect to the Eu- 
clidean time is imposed on the Bose (Fermion) fields. 
This results in the discrete Matsubara frequency and the 
summation over it is involved. In the real-time formal- 
ism, on the other, we deal with the real (Minkowski) 
time directly, at the cost of extending the theory to the 
complex time. The complex time path has some arbi- 
trariness under some constraint, and one possible choice, 
which is usually employed, is the one composed of the 
three paths: G+ = [— oo : oo], G_ = [oo — it] : — oo — iry] 
and C vcr ticai = [— oo — irj : — oo — i/3], among which the last 
piece eventually decouples from the theory. Correspond- 
ing to the two time paths of G+ and G_, we have two 
kinds of interaction vertex, and the propagator becomes 
a two-by-two matrix form. 

CGC is formulated in the light-front field theory [I?], 
HH, so that we need the light-front field theory at finite 
temperature when considering the thermal CGC. The 
light-front thermal field theory was initiated in Ref. [4£| 
and formulated in Refs. [U M, HI, [H, H, In the light- 



front thermal field theory, the (imaginary and real) time 
variable now refers to the light-cone time x + . It is noted 
that the partition function is given by Tr[exp{ — (P + + 
P~)/\/2T}] rather than Tr[exp{-P"/T}] where P~ is 
the light-front Hamiltonian. Also, the theory has non- 
trivial dependence on the frame we choose. It is argued 
that it is impossible to formulate the theory in the frame 
where the heat bath moves at the speed of light. The the- 
ory formulated in the rest frame of the heat bath turns 
out to be the simplest. In the following, we will work in 
the rest frame of the heat bath, which is also the phys- 
ically relevant frame for describing an actual heavy ion 
collision that we are interested in. 

For our purpose, we employ the real-time formalism 
rather than that based on the imaginary-time. We have 
two reasons: (1) As is argued in Ref. [35|, CGC is formu- 
lated properly on the complex time contour C — G+UG_ . 
In the real-time formalism, it is possible to obtain the fi- 
nite temperature theory just by adding the path C ve rticai 
to G. (2) When we derive the RG equation which gives us 
the BFKL and JIMWLK equations, we have to perform 
the shell-momentum integration over the energy variable 
p~ . In the imaginary-time formalism, we have to do the 
analytic continuation from the discrete Matsubara fre- 
quency to the continuous energy variable, while in the 
real-time formalism, we can treat the continuous p~ di- 
rectly. 

The bare real-time thermal gluon propagator in the 
light-front theory G^ a ^ was given in Ref. [H^], where 
fi, v = ±, i (i = 1, 2) are the Lorentz indices and a, (3 — 
± are the indices associated with the contours G+ and 
G_, and the color indices are suppressed. We put the 
indices a, (3 in parentheses in order to distinguish them 
from the Lorentz indices which also take the values of ±. 
In Appendix [XJ we explicitly write down the non-zero 
components of G^ a m in the light-cone gauge (A+ = 0), 
which will be used in the following discussions. 

As is emphasized in Ref. [3o| , it is important to specify 
the prescription for the pole in l/p + . If we use an inap- 
propriate prescription, we are in a trouble as the soft 
gluon emission becomes non-local: There appear smaller 
x gluons in negative x~ even if the color source lies in 
positive x~ . This has a serious consequence especially 
for deriving the JIMWLK equation. We use the retarded 
prescription as in Ref. [351 ] , in which all the radiated glu- 
ons come in positive x~ definitely as it should be. 

It is noted that there is still an arbitrariness for the 

pole prescription for the G^ ) and G( |_) components 

even if G( ++ )(= G?_0 is specified for which we use 
the retarded prescription. In this paper, we take the 

same prescription for the G(_| ) and G( |_) components 

as for G( ++ ), as is shown in Appendix lAl In fact, Ref. 
[55l | proposes a more symmetric prescription in which, for 
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instance, the Gw_| ^ component has the form of 
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In the symmetric prescription, the pole prescription is 
replaced by its Hermite conjugate when it is multiplied 
by Gq. It can be shown explicitly that in the symmet- 
ric prescription, there arises a contribution of the gluon 
distribution in negative x~ , while in our prescription, 
the radiated gluons appear only in positive x~ . Thus 
our prescription is the appropriate one for calculations of 
CGC. 



III. COLOR GLASS CONDENSATE AT FINITE 
TEMPERATURE 

In this section, we discuss how the CGC formalism 
is extended to finite temperature. Let us begin with 
recapitulating the theory at zero temperature [35l |56| . 
CGC is the high gluon density matter that appears in 
the high energy limit of hadrons and is relevant to the 
hadron scatterings in the Regge kinematics. Specifi- 
cally, we can imagine the high energy hadron-hadron 
(or nucleus- nucleus) collision or Deep Inelastic Scattering 
(DIS) where a high energy virtual photon collides with 
a hadron. The generating (or partition) functional for 
CGC is given by 

Z\j} = J VpW A \p]Zl l \p] I B^KV M "' /i,Jl , 

which contains the generating functional for the soft 
gluon field at a fixed color charge distribution p; 

Z A [ P ,j] = Zl l [p] [ A VA%5(Ai) e^M/i^ (2) 



with Za [p] = Z\ [p, j = 0] . The frozen color charge p 
simulates the effect of the fast partons with momenta 
|p + | > A+, and the charge distribution is specified by 
the weight function Wa[p] into which all the dynamics 
of the fast modes is encoded. Let us recall how we can 
calculate the gluon distribution function within the for- 
malism. There are two ingredients: 

(i) We solve the classical Yang-Mills equation to some 
specific charge distribution, assuming the presence 
of the strong gluon field at the small x region of 
interest. 



(ii) The gluon distribution function is obtained by av- 
eraging the classical solution A l (x = (x~ ,x±))[p] 
we have obtained over p with the weight function 
W A [p\. 

We note that the classical solution „4 J (:?)[/>] we need 
is a static, time-independent one, which is necessary for 
computing the gluon distribution function. 

If we lower the cutoff A + by integrating out the semi- 
fast gluons a M (6A+ < \p + \ < A + ), their dynamics is 
rcnormalized into Wi,a[p] to give the RG equation for 
Wa[p]- Convoluted with the classical solution, this RG 
equation for W\ [p] finally gives us the evolution equation 
for the gluon distribution function, that is, the BFKL and 
JIMWLK equations. The RG equation for Wa[p] reads 



dW T [p] 

dr 



2 8p T {x)8p T {y) 



[WrXxy 



Spr(x) 



[W T a x 



where 



a a (x) = (5p a (x)), 
Xab(x,y) = (5p a (x)6p b (y)}, 



(3) 



(4) 
(5) 



and r = ln(P + /A + ) = ln(l/a;) is the rapidity variable. 
a(x) and y) are the one- and two-point functions, re- 
spectively, of the induced color charge which is generated 
by the semi-fast gluon fluctuation and is defined by 



SJ a Jx) 



ss 



6AZ(x) 



A+a 



5 2 S 



6A£(x)5AZ(y) 
1 S 3 S 



<{y) 



2 6M(x)8A%(y)5A${z) 



al{y)a x c {z) (6) 



with fx = — . The derivation of the RG equation is re- 
duced to computation of <x(x) and x(a:,y). 

Now we consider the extension to finite temperature of 
the formalism. First, we note that the classical solution 
at finite temperature is identical with that at zero tem- 
perature: The thermal effect does not alter the classical 
solution within the classical approximation. This can be 
seen as follows. The thermal effect in finite temperature 
field theory is brought about by the modification of the 
time contour on which the theory lives. As we noticed, 
the classical solution needed in the CGC formalism is the 
time-independent one. Thus, the classical solution at fi- 
nite temperature is provided by the same one as at zero 
temperature. 

All the thermal fluctuations are contained in the cor- 
relators ct(x) and x(a;,y), the evaluation of which is all 
we have to do. The action in the real time formalism at 
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finite temperature is given by 

S[A,p] = Sym + Sw, (7) 

Sym = ~ [ d^x-F^F^, (8) 

S w [A-,p] = JL J d 3 fTr [p(x)W c [A-]{x)] (9) 

with the time-ordered Wilson line 

W c (x) =T c exp^ig J^dzA-(z,x)X , (10) 



where C denotes the complex time contour given by C = 
C + U C_ U C vcrt icai, f° r which the last piece is eventually 
decoupled from the theory (45[. With this action, we can 
calculate the induced charge 5p defined in Eq. ©, which 
reads 



5p a (x)=6p2Hx) + 6pW(x), (11) 



where 



2igT:;{x)a lc {x) 



gPac{x) 



dy+ ) (x+|PV^|y+)a c -(y+ ) ,f) + i 



dy^a {y^yx) 



(12) 



5p^\x) = 9 2 f abc [d + at(x)] aftx) 



dz+ +) a c (y+ +) ,x)a d (z+ +) ,x) 



{0(x+ y+ +) )e(y+ +) z+ +) )Tr{T a T c T d T b ) + 9(y+ +) - z+ +) )6(z+ +) x + )Tr(T a T b T c T b ) 
+0{y+ +) - x+)9(x+ - z+ +) )Tr{T a T d T b T c )} 

pOO 

dy+_ } I dz+ +) a c -(y+_ y x)a d -(z+ +) ,x)6(z+ +) - x+)Tr(T a T b T c T d ) 
dz+_ } a c -(y+ j, x)a d -(z+_ >, x)6(z+_ ) - y+_ } )Tv(T a T b T c T d ) 



— OO 

oo 



/OO p CO 

dy+) / dz+ ) a c -(y+ ) ,f) a d -(z+ ) ,x)0(x+-z+ ) )Tr(r a T' i r b T c ) 



(13) 



It is noted that the terms that involve the integration over 
the time variable on C_ (i.e., yf^ and ztO represent the 
thermal effects, while the other terms are already present 
at zero temperature. 

Having obtained the charge fluctuation, we can now 
express a(x) and x(x,y) in terms of the full Feynman 
propagator G^^(x, y) of the semi-fast gluons. To the 
lowest order of a s , we have 



where 



at{x) 



g 2 ir bc d+Gt ++)cb (x + ,x;y + ,y) 



x=y 



-g 2 d+Tr T a Gf ++) (x+,x;y+,y) 



(15) 



a a {x) = a a l {x)+& a 2 {x) 1 



(14) 
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+ 



dy (+) / dz (+) G (++)d C ( 2: (+)' :z; ;y(+)' a; ) 

) J — OO 

{9(x+ - y^ +) )e(y+ +) - z+ +) )Tr(T a T c T d T b ) + 9(y+ +) - z+ +) )6(z+ +) - x+)Tr(T 
+%+ +) ~ x+)6(x + - z+ +) )Ti(T a T d T b T c )} 

dz+ +) G(-_ )dc (z+ +) ,x'; y+_ v x)6(^+) ~ x + )Tr(T a T b T c T d ) 
dz+_ ) G ( -Jl )dc (z+ j, x; vt_ v £)B{*U ~ V U )Tr(T Q T b T c r 
<! V , / &z+ +) Gj-_ )dc {z+ +) ,x-y+_ ) ,x)6{x+ - z+ +) )Tr(T a T d T b T c ) 

I 



(16) 



and 



where 



x(%, y) = Xi(x, y) + x%(x, y), 



(17) 



Xi{x,y) = 4ig 2 ^ +i (x)G^ ++) (x,y)^(y), (18) 



lx2(^y)=^ +l (.x)|-2 J dw+ +) G\l +) (x+,x-,w+ +v ff){w + \PV^:\y + )+i J dw+_ ) G\-_ ) (x+,x; w+_ p{g) 
+ p(x)! [ 2 J dz+ +) (x+\PV^\z + )Gll +) (z+ +) ,x;y+,y) + i J dz+_ } G^ +) (z+ } , x ; y+ , £)} ^ (y) 



+ p(x)ii J dz+ +) J Aw+ +) (x+\P\— \z+)G-- +) (z+ +) ,x;w+ +r y)(w + \PY— \y+) 



+ 2j dz (+) J dw t-)( x |PV laf |z ) G (+-)( Z (+)'^ W (-)^) 



r d. 



z (-) J dw (+) G (-+)( z (-)'^ w (+)'^)( w l pv ^|y ) 



(19) 



i 

The diagrammatic representations are the same as derive the thermal BFKL equation within the finite tem- 
given in Ref. [IH, although the internal time variables perature CGC discussed in the previous section, using 
can now take values on C_ , giving rise to the additional the same approximation. All we have to do is to evalu- 



terms corresponding to the thermal effects. 



IV. BFKL EQUATION AT FINITE 
TEMPERATURE 



ate (7 and x m the approximation, for which we need the 
free propagator G^Y^s presented in Appendix [A] The 
discussion in this section is proceeded in parallel with 
that given in Sec. 5 of Ref. [35j |. 



It has been shown that the CGC formalism at zero 
temperature reproduces the BFKL equation in the weak 
source approximation (3l| . In the present section, we 
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A. Evaluation of a in the weak source limit 

We first consider o\ . Because of the color trace in Eq. 
(fT5)l . the vacuum contribution disappears. The contri- 
bution linear in p that has logarithmic enhancement is 
given by a direct insertion of p, which reads 



a i (x) = — g 2 / d 4 z / d 4 u 
Jc Jc 

xG l -(x - z)Tl c \z, u)8+Go l (u - y)\ 

where 



S 2 S W 



8A-{z)5A-{u) 



~PabS(z- U) {9 C {Z^ 



(20) 



u+) - 6 c (u+ - z+)} 
(21) 



after the rescaling dp — > ^i5p. We note that because the 
time contour C has the two pieces, there arise four terms 
for a x in Eq. (JUJ. 



In the following, we directly calculate the two- 
dimensional charge density cr a (x^) = / dx~<T a (a;), which 
is actually what we need in the RG equation. In momen- 
tum space, Eq. (|20|) is written as (with q=p — k) 



d-i (£) =- o.9 2 



dp 
"2V 



d 3 p 



d 3 fc 



(2tt) 3 y (2tt) 



p(g)<j-2G^ ++) (p-,p) 



PV- 



fc+G 



o(++) 



CP - ,*:) 



+ 2 7 rW(p-)G^ ++) (p-,p)fc+G 7_ +) (p-, k) - 2iti6(p-)&- {+ _ ) (p-,0)k + G-{ ++) (p-,k) 



+2G i 7 + _ ) (p- J ^) 



PV — 



" + G ! o(-+)0 J 



(22) 



As in the case of zero temperature, we can find that k + 
could be replaced by p + or (k + + p + )/2, by noting the 
symmetry of the integrand under the simultaneous ex- 
change k <-> — p and p~ «-> — p~ . If we take the strip 
restriction on p~ (A - < \p~\ < A~/b), the second and 
third terms in the curly brackets in Eq. (|22p disappear 
due to the delta- function of p~ . After a straightforward 
calculation similar to that at zero temperature, we obtain 



a<!(x x ) = 
d 2 p A 



9 2 N C 
2tt 



A~/6 



A- 
X 



dp-i(l + 2n B (^ 



PA 



fcj 



(2tt) 2 y (2tt) 



pikj 



(23) 



where p a (q±_) = p a (q + = 0,q±) and nB(x) = 
l/(exp(/3x) — 1) is the Bose distribution. It should be 
mentioned that the strip restriction on p + (bA + < \p + \ < 
A + ) leads us to the same result as for the p~ restriction, 
meaning the implicit symmetry under p + <-> p~ of the 
integrand in Eq. (f2"2")l . Fourier transforming Eq. (|2"3"]) . we 
finally have 



al{q±) = 2N c a s 



A- /b 



d,--L(l + 2, B (^ 



x P a (q±) 



d 2 P± 



<ii 



{2nY\2p\{p 1 _-q 1 _Y p\ 



We note that at zero temperature, the integral over p~ 
gives rise to the logarithmic enhancement factor ln(l/6), 
while in the present case at finite temperature, the inte- 
gral cannot be performed analytically. 

One comment is in order. We could calculate a(x) in 
the x~ representation, explicitly showing the longitudinal 
structure. As we have noted, it is possible to verify that 
a(x) has contributions only at positive x~ in the non- 
symmetric prescription presented in Appendix fAl while if 
we use the symmetric prescription mentioned in SecHH it 
turns out that there appear contributions at both positive 
and negative x~ . Thus our non-symmetric prescription 
is valid for calculations of CGC. 

Similarly, we can compute erf (q±) to be 



a$(q ± ) =2N c a s J^ ' dp-^ (l + 2n B (jL 

x/(9i)f^T. (25) 



(2tt) 2 P 2 ± 



which cancels the tadpole term \jp\ in ax- Combining 
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the two pieces, we finally obtain 



=N r a R 



\~/b 



d %-( 1 + 2nB (S 



X p a (q ± ) 



d 2 p± 



(27r) 2 p2>x-9±) 2 



with 



\-/b 



X^(x,ff)=4a s J^ dp- ±. (l + 2n B 

r|2 p ipx-(x±-yx) 

i ^ * {p^)p(f) + (J*W) 



(2tt)2 

(26) -i (J^W) . 



(31) 



B. Evaluation of x in the weak source limit 

We start with \i given in Eq. (fl"8|) . which, in lowest 
order in p, reads 

xf(x,y)=4ig 2 ^(x)G^ (++)cd (x,y)^'(y) . 

(27) 

As in the zero temperature case, we can manipulate 

G o\++) ( x > y) t0 be 



u o(++ 



(2tt)2 




iPx-Ox-I/x) 



where we have imposed the strip restriction on the p + 
integration, though it is possible to show that the same 
result can be obtained when we take the restriction on 
p~. Up to the leading logarithmic accuracy (LLA), we 
can replace p\ in the Bose distribution by Q\ which is 
some typical transverse momentum. Then the factor in 
the parentheses can be factorized out of the p± integral. 
Furthermore, as we will see later, we assume the temper- 
ature is comparable with P~ jx — p~ with P~ being the 
nucleon light-cone energy, which is much larger than p + 
so that the term p + in the argument of the Bose distri- 
bution can be neglected. It is noted that the p + term is 
also much less than the term Q 2 l /2p + « p~ in the argu- 
ment for the on-shell excitation. After a simple change 
of variable, we obtain 



G 



o(++)( x >y) — 



IS 11 \ 
Ja- 



2?r 2p~ 



2nB (7! 



y±) 



(29) 



which gives us 



x f(x,y)=2g 2 



\-/b 



2vp-\ B \V2 



x 5™(x ± - yi _)T^{x)Tt b \y). (30) 
In a similar way, X2 b can be calculated, providing us 



C. Thermal BFKL equation 

Now we derive the thermal BFKL equation, using the 
RG equation ([3]). As we have mentioned, contrary to 
zero temperature, the enhancement factor ln(l/6) is not 
manifest at finite temperature, for which we have to be 
careful. If we make the factor explicit, Eq. ((3| is ex- 
pressed as 

^M = Hm^ r 
dr b^i In i 



1 



2 Sp T (x)Sp T (y) 



[WrXxy] 



8p T (x) 



[W T a x _ 



(32) 



We are interested in the evolution of the diagonal element 
of the two-point function (pp) T , 

( Pa (kx)Pa(-kx))r « *? ± < (k ± ) | ) T = v(x,k 2 ± ), (33) 

which is called the unintegrated gluon distribution. 
By multiplying Eq. (|32|) by J d 3 xd 3 yexp{—ik± ■ (x± — 
yx)}p(x)p(y), and by using Eqs. ([26]) and §1} for cr (0) 
and ) we obtain 



0_ 

dr 



ip(x, k 



lim — -r 

6-1 In \ 



a s N c . 



hr lb 



dp 



A- 



1 + 2n B 



V2 



dV" 



T<p(x,p 2 j_) 



P±(P± - k±) 2 

\S d2 ^Piipl-k^ kl 



(34) 



We now discuss how the logarithmic enhancement fac- 
tor can appear from the prefactor 



A~/6 



dp 



1 



1, 



in Eq. (f3~4"|) . We note that in the present framework, 
the temperature T is entirely a free parameter which can 
take any value. It is useful to divide three temperature 
regimes as follows: (i) T <C A~/b, (ii) T — A~/b and (hi) 
T ^> A~ /b. It is noted that as is indicated in Eq. (j3~4"| . 
we take the limit of b — > 1 in the end so that the three 
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regimes will be translated to (i) T <C A - , (ii) T ~ A - and 
(iii) T ^> A - after all the calculation. In the regime (i), 
the thermal distribution function is exponentially small 
in the integral shell so that the thermal effect plays no 
role up to LLA and we find that the BFKL equation re- 
duces to that at zero temperature. In the regime (ii), the 
distribution function is substantial at the upper bound- 
ary. The thermal effect contributes to the logarithmic 
enhancement. This is so because the logarithmic factor 
comes from the upper boundary of the integral. In the 
regime (iii), the temperature is much higher than the 
shell energy. In this case, we find that there appears 
a power enhancement rather than the logarithmic one. 
This means that the BFKL equation which resums the 
leading logarithmic terms makes no sense. Physically, 
the eikonal approximation breaks down in this regime. 

In the following, we restrict ourselves to the regime (ii). 
Although the evolution in the regime (iii) is interesting, 
it is beyond the scope of this work. 

Let us now extract the logarithmic factor. For that 
purpose, we need to be careful. As is seen in Eq. (j34| . 
we eventually take the limit of b — * 1 that makes the in- 
tegral shell infinitesimally small, leading us to the BFKL 
equation being a differential equation. Thus, we are es- 
sentially considering the temperature close to A~ /b with 
6=1, that is, T w A - . However, in order to extract 
the logarithmic factor, we have to temporarily take the 
limit of b — > before taking b — > 1, as in the case of the 
vacuum BFKL equation [351 ] . This means that the upper 
boundary of the integral virtually moves toward infinity. 
Thus, in this intermediate stage, it is necessary to impose 
the constraint of T as A~/6 in order to keep the temper- 
ature staying at the upper boundary, or in order to keep 
ourselves staying in the regime (ii). With this constraint, 
we expand the prefactor as 




logarithmic factor. It is easy to take the operation of 
lim(,_>i 1/ ln(l/6), finally leading us to the thermal BFKL 
equation 



ip(x,k\) 



a s N c 



X { I d V n 



6% /2Tx — 1/ 



= a_i In - + a + a\b + a 2 b 2 + ■ ■ ■ (35) 
b 



-\I^ A^-^) M x M- (37) 

where we have used the relation A~ = P~ jx with P~ 
being the nucleon light-cone energy. The equation is 
quite simple in the sense that the only difference from 
the BFKL equation at zero temperature is the Bose en- 
hancement factor in the prefactor. This can be given a 
natural physical interpretation that the emission of the 
softer gluons from the harder ones is promoted because 
of the presence of the thermal gluons all around. We 
conclude that at finite temperature, the growth of small 
x gluons can be more rapid and the saturation regime 
can be reached sooner than at zero temperature, if the 
temperature is so high as to be in the regime (ii) , namely 
T~P~/x. 

Finally, we comment on a different thermal BFKL 
equation that has been derived in Ref. 57]. In their 
equation, the thermal modification appears in the kernel 
in transverse coordinates, while our equation involves the 
overall Bose prefactor. The difference is due to the fact 
that they assume the thermalization for modes in the t- 
channel at the center-mass frame, while in our framework 
it is for modes in the s-channel. Their thermal BFKL 
equation would be useful for investigating Pomerons at 
finite temperature as argued in Ref. [5jJ , and ours would 
be relevant to the initial state of heavy ion collisions and 
the jet quenching. 



SUMMARY 



with b -C 1. In order to find the coefficient a_i, we 
differentiate the both side with respect to b to have 



1 



e V2bT — i 



= a_i — a\b — 2a 2 b 2 + 



If we take 6^0 with keeping A jbT 1, we obtain 



a_i = 1 + 



(36) 



g \/2bT — \ 



Thus, the thermal effect can contribute to the logarithmic 
enhancement in the regime (ii). 

Now that we have extracted the logarithmic enhance- 
ment factor, we turn to take the limit of b — > 1. We note 
that in this final stage, we should abandon the constraint 
of T w A _ /6, the purpose of which is just to extract the 



We have derived the thermal BFKL equation within 
the formalism of Color Glass Condensate (CGC). CGC 
at finite temperature would be relevant to the initial state 
of ultra-relativistic heavy ion collisions and to understand 
the anomalous cross-section of energetic partons passing 
through the quark-gluon plasma (QGP). We have dis- 
cussed the extension of CGC to finite temperature and 
found that the classical solution is not modified at fi- 
nite temperature and that all thermal effects are included 
in the induced charge correlators &{x) and x(£,y) ap- 
pearing in the renormalization group equation for the 
weight function Wa[p], In the weak source approxima- 
tion, we have derived the thermal BFKL equation by 
explicitly evaluating the charge correlators in the real- 
time formalism. We have specified the pole prescription 
for l/p + in the free gluon propagator, as is presented in 
Appendix [XJ We have employed the retarded prescrip- 



tion for the and G( )(= components, as 

well as the non-symmetric prescription for the G( H ) and 

G(_ + ) components. The prescription is justified by the 
fact that the gluon radiations occur definitely at positive 
x~ . 

We note that there are the three temperature regimes: 
(i) T < P~/x, (ii) T ~ P~/x and (iii) T > P~/x. 
In the regime (i), thermal effects disappear up to LLA 
and the evolution equation reduces to the vacuum BFKL 
equation. In the regime (iii), the temperature is too high. 
We have a power enhancement rather than the logarith- 
mic one so that the evolution is described by an equation 
different from the BFKL-type equation which resums the 
leading logarithmic terms. This regime is beyond the 
scope of the present paper. In the regime (ii), the tem- 
perature matches the light-cone energy of soft gluons we 
are thinking of, giving rise to the logarithmic enhance- 
ment. Thus, in this regime, it is meaningful to consider 
the thermal BFKL equation. The resulting BFKL equa- 
tion indicates that the thermal effect shows up as the 
Bose enhancement of the soft gluon emissions due to the 
surrounding thermal gluons. We have concluded that 
at finite temperature, the growth of small x gluons can 
be more rapid and the saturation regime can be reached 
sooner than at zero temperature. 

To complete, the full non-linear JIMWLK equation at 
finite temperature will be derived in the future. 



iG o(++)(P) = 



iG, 



o(+-) 



(P) 



p+ - 


ic 


i 




P 




p+- 


ic 






P+- 


ic 







[iG Q (p) + 27rn B (\p°\)6(p 2 )], 
2tt [0(-p°)+n B (\p°\)]6(p 2 ) 



iG-^ip) = * 2* [0 (pP) + n B (|pP|)] S (p 2 ) 



p + + ie 



iG 



o(++) 



(A2) 



(p) = S i i[iG (p) + 27rn B (\p°\)S(p 2 )], 
iGg + _ } (p) = Mir [0 (-p°) + n B (\p°\)] S (p 2 ) , 
iGl_ +) {p) = Mir [9 (p°) + n B (\p°\)] 6 (p 2 ) , 
K\-)(P) = 511 ["iGS(p) + 27m B (\p°\) 5 (p 2 )] , (A3) 
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APPENDIX A: REAL-TIME THERMAL GLUON 
PROPAGATOR IN LIGHT-FRONT FIELD 
THEORY 

In this appendix, we write down explicitly the real- 
time thermal gluon propagator in light-front field theory 
with the light-cone gauge (A+ = 0). Note the retarded 
prescription for G l ~ and G~ l components. 



iG -- +) (p) = PV^- [iGo(p) + 2nn B (|p°|) S (p 2 )] , 
= Pv| r 2 7 r [6 (-/) + n B (|pP|)] 8 (p 2 ) , 

iG - ( - + )W = PV 7T 2 - i° (P°) + (\P°\)] 5 (V 2 ) , 
iG-^ip) = PV^- [-iGS(p) + 2nn B (\p°\) S (p 2 )] , 



(A4) 



where 



iG o(++)W 



p + + if 



p + + ie 



[iGo(p) + 27Tn B (\p°\)S(p 2 )] 
2tt [e(-p a )+n B (\p a \)]S(p 2 \ 



i^_ +) (p) = ^-2. [0 (pP) + n B (\p°\)} S (p 2 ) , 

iG^__ } (p) = [-iGS(p) + 2nn B (|p°|) S (p 2 )] , 

(Al) 



G Q (P) 
n B (x) 



1 



p A + ie ' 



1 



2 \p+ — ie p+ + ie 
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